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ABSTRACT 

Input/Output Cost Analysis involves decompositions of the quadratic cost function 
into contributions from each stochastic input and each weighted output. In the past, these 
suboptimal cost decomposition methods of sensor and actuator selection (SAS) have 
been used to locate perfect (infinite bandwidth) sensors and actuators on large scale 
systems. This paper extends these ideas to the more practical case of imperfect actuators 
and sensors with dynamics of their own. NASA’s SCOLE examples demonstrate that 
sensor and actuator dynamics affect the optimal selection and placement of sensors and 
actuators. 
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1.0 INTRODUCTION 


The objective of this paper is to develop and evaluate a method for the selection of 
sensors and actuators in the control of finite-dimensional linear systems using imperfect 
sensors and actuators — devices which do not provide instantaneous responses, but have 
nontrivial dynamics of their own. In addition, the plant noise and the measurement noise 
is assumed correlated. This important case allows the use of accelerometers as sensors 
(this always yields correlated plant and measurement noise). Application of the 
generalized method to practical control problems demonstrates that correlatedness of the 
noise and the dynamics of the actuator and sensor devices can significantly affect the 
optimal selection of both the number and location of sensors and actuators. 

Consider as a starting point the following familiar dynamic system model: 


X P = A p x P (t) + B p [f(t) + w(t)] 

(1.1a) 

y p (0 = C p x p (t) , z(t) = M p x p (t) + v(t) 

G.lb) 


E { w(t)w T (x) ) = 8(t-t)W , E {v(t)v T (t)} = 5(t-x)V , E{ w(t)v T (t)} = 5(t-x)U(l.lc) 

where x p e R *, fe R"\ w e R w , z,v, € R" 1 and (A p ,C p ) observable, (A p ,B p ) 
controllable and (A p ,Mp) detectable. The vectors w(t) and v(t) are respectively zero 
mean white noise characterizations of the actuator and sensor noise. 

In control of large space structures, the locations of sensors and actuators becomes a 
critically significant "degree of freedom" in control design [14, 20]. Among over 60 
more recent contributions to the SAS problem, only [4], [7], [10], [11], and [12] consider 
noisy actuators (W, V nonzero). In all cases, the disturbances are modelled as Gaussian, 
white, and uncorrelated (W, V diagonal, U = 0). Most of the SAS literature takes no 
account of actuator or sensor dynamics. Two exceptions are McClamrock [19], and 
Howell and Baxter, [6]. In [1] the authors extend the cost decomposition approach [2] to 
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accommodate noise correlation between sensor and actuator noise sources (W, V not 
diagonal, U * 0). A key conclusion in [1] is that the proper sensor/actuator selection and 
placement can be drastically affected by noise correlation. For example, the deletion of a 
noise source (by making an actuator or sensor noise free) may degrade performance 
contrary to the usual expectations when noise sources are uncorrelated. 

Very fast actuator dynamics may be neglected in stability considerations, [9]. A 
more thorough discussion of the effect of actuator dynamics is given by Goh and 
Caughey [8]. The analysis of [8] and [9] demonstrates that plant frequencies occurring 
above the actuator bandwidth can lead to closed loop instability, even for co-located 
sensors and actuators. Goh and Caughey do not address the problem of selection of 
dynamic actuators. That is the goal of this paper. 

This paper is organized as follows. First the system model is augmented to include 
sensor and actuator dynamics. The closed-loop input and output costs are then developed 
for the fully augmented system, and they are used to define expressions which reflect (he 
effectiveness of each dynamic actuator or sensor in minimizing the cost function. 
Finally, the method is illustrated by application both to small scale numerical examples 
and to NASA’s SCOLE flexible space structure model. It is found that in the selection of 
noisy actuators and sensors, finite dynamics can significantly affect selection results. 

2.0 MODELING DYNAMIC ACTUATORS AND SENSORS 

In [2] the results of Closed-Loop Input/Output Cost Analysis (CIOCA) were 
developed and applied to the Sensor and Actuator Selection problem (SAS) for systems 
of the form (1.1) under closed-loop control. In [1] the control f(t) is the vector of optimal 
state estimate feedback controls: 
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(2.1a) 


f(t) = Gx c (t), G = -R _1 BjK. 

X c = Vc(‘) + B p f W + F t z (0 - MpX C (t)l • F = [ pM 7 + BpUJV- 1 , (2.1b) 

0 = KAp + A p T K - KB p R- 1 B p t K + C p T QC p (2.2a) 

0 = [A p - BpUV-'MpJP + P[A p - B p UV _1 Mp] T - PMjJv^MpP (2.2b) 

+ B p WB p T -B p UV- 1 U T Bp T 
which minimizes the cost function 

V =£„{Hy p (t)ll 2 Q + llu(t)ll 2 R} , E^lim E[-J (2.3) 

t— 


where x c e R n " is the vector of state estimates. The conclusion from [1] for this problem 

(1.1), (2.1) (2.2) is that when U*0, the sensor/actuator selection results can be 
drastically different. 


2. 1 Adding Actuator Dynamics 

First the system (1.1) is augmented to include stable, observable, controllable 
actuator dynamics of arbitrary order. 

x a = A a x a f B a (u+w u ) , (A a ,B a ) controllable (2.4a) 

f = C a x a , (Ag.CJ observable, x a e R"*, fg R"* (2.4b) 

Figure 2.1 presents schematic representations for actuator models of varying degrees of 
complexity; Figure 2.1a represents the non-dynamic actuator, while Figure 2.1b 
represents the general model for a dynamic actuator with white noise. Note that for the 
non-dynamic actuator the noise w(t) is purely additive with the input u(t). In the case of 
dynamic actuators the analyst may consider the actuator’s output (into the system) to 
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Figure 2.1a 



Figure 2.1b 


Figure 2.1: Actuator Models 
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include additive actuator output noise Wf(t), or actuator command noise w u (t) which is 
filtered by the dynamics of the actuator, or both. Both types of noise are assumed 
possible in our development. 

Augmenting the system states x p of the original system (2.1) with a vector of 
actuator states x a , we obtain: 

x = Ax + Bu + Dw , y = Cx , z = Mx + v = z p (2.5a) 
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where f = C a x a , (A a , C a ) is observable, and Re 
is controllable. 




<0, i=l,2, ••• n a . (A a ,B a ) 


First note that since x p is observable from y p , (i.e., (A p ,Cp) is observable) and x a is 
observable from f (i.e., (A^CJ is observable) then from the definitions (2.5) x must be 
observable from y, that is: 


(A,C) is observable 


(2.6a) 


Also note that the actuator dynamics are assumed stable, so that the system (2.1) has 
not been augmented w'ith any unstable states. Therefore, detectability of (A p ,M p ) 
together with stable A a yields 

(A,M) detectable. (2.6b) 
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Finally, Theorem 1 states the conditions for controllability of the system (2.1) 

augmented with actuator dynamics (2.5). Proof of the theorem is contained in the 

Appendix. 

Theorem 1 

Consider the controllable system 

x p = ApXp + B p (f+w f ) , (A p ,Bp) controllable (2.7a) 

x p 6 R™* (2.7b) 

augmented with controllable and observable actuator dynamics of arbitrary order 

x a = A a x a + B a (u+w u ) , (A a ,B a ) controllable (2.7c) 

f=C a x a , (A a ,C a ) observable, x a eR nx *, fe R nu (2.7d) 

to form the composite system 

x = Ax + Dw + Bu fl.lt) 



'A p B p C a 


'B p O' 


0 


X P 


Wf 

A = 

.0 A a 

, D = 

0 B a 

, B = 

Ba 

, x = 

/a 

, W = 



The system siates x p are controllable from u(t) if the number of poles minus the 
number of zeros is the same for each individual actuator’s transfer function. 

Remark 1: Note that full controllability of the augmented-system state vector x is 

not guaranteed under the conditions of the theorem. 

Remark 2: The conditions of the theorem are always met for first order dynamic 

actuators, (assuming no direct input/output "feedthrough" for the 
actuators), since each actuator will have one pole and no zeros. 
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Remark 3: 


As long as the original system states Xp are controllable through some 
minimum set of actuators meeting the criteria of the theorem above, then 
controllability of Xp will be maintained with the addition of actuators of 
any order and any number of transfer zeros. 

Remark 4: Finally, note that the usefulness of the theorem stems from the fact that 

by meeting certain mildly restrictive conditions, the actuator dynamics 
can be guaranteed not to destroy controllability of the original system 
states, regardless of the pole/zero locations of the plant. 

2.2 Adding Sensor Dynamics 

Next the system (2.5) is augmented to include stable, observable, controllable 
sensor dynamics of arbitrary order. 

x s = A s x s + BjOlpXp+v^), (A S ,B S ) controllable (2.8a) 

z = C s x s + v out , (A S ,C S ) is observable, x s g R 11 * (2.8b) 

Re[X i(A s )]<0, i = { 1, 2, • • • n s ) (2.8c) 

Figure 2.2 presents schematic representations for actuator models of varying degrees of 
complexity; Figure 2.2a represents the non-dynamic sensor, while Figure 2.2b represents 
the general model for a dynamic sensor with white noise. Note from both eqn (2.8) and 
Figure 2.2b that (similarly to the case of actuator dynamics), adding sensor dynamics 
leads to the possibility of both sensor input noise and sensor output noise. The i^ 
sensor’s input noise (v in )j is filtered by the dynamics of the i^ sensor, while the output 
noise (v^J, is purely additive with the sensor output. Both types of noise are assumed 
possible in our development. 
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Figure 2.2b 


Figure 2.2: Sensor Models 
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The fully augmented system equations have the following form: 


x = Ax + Bu + Dw 
y =Cx 
z = Mx + v 


(2.9a) 

(2.9b) 

(2.9c) 


* T = [x p T , x a T , x s T ] , y T = [y p T , f T l , w T = [w f T , w u T , vj] , v = v, 
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, M =[0 0 C s ] 


The response y p (s) of the plant to the input f(s) is given by 

y p (s) = Hp(s)f(s) (2.10a) 

where 

Hp(s) = Cp(sI-A p ) _1 B p . (2.10b) 

is the plant transfer function. The response Mx(s) of the actuator/plant system (2.9) to 
the input u(s) is given by 

Mx(s) = H(s)u(s) , (2.11 a) 

where 

H(s) = M(sI-A) -1 B. (2.11b) 

Finally, the response z(s) of the sensors to an input Mx(s) is given by 
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z(s) = Hs(s)Mx(s) , 


(2.12a) 


where 

H s (s) = C s (sI-A s r 1 B s (2. 1 2b) 

is the transfer function for the sensor dynamics. Minimal systems are controllable and 
observable. Thus, given minimality of the plant/actuator system [(A,B) controllable and 
(A,M) observable], then measurability of the full augmented system is guaranteed 
[(A ,M) observable] if there are no pole/zero cancellations between H(s) and H s (s). 

2.3 Defining the Cost Function 

With the properties of the augmented system established, optimal control design for 
the augmented system is now considered. Recall that the standard LQG cost function 
(2.3) for the unaugmented system (2.1) includes a penalty on the output regulation error 
y(0, as well as a penalty on the control energy u(t). However, in the augmented system 
(2.5), while the actuator command is given by u(t), the actuator response f(t) (contained 
in the augmented output vector y) is distinct from u(t) due to actuator dynamics. A true 
measure of control energy is more appropriately stated in terms of a weighted sum of the 
variances of fit) rather than of u(t). It can readily be shown, however, that even in the 
presence of a weighting on the actuator outputs, fit), some nonzero weighting on the 
actuator inputs u(t) is necessary to avoid an infinite gain solution to the optimization 
problem. For this reason, and in view of the relation of fit) to the design goals as 
discussed above, minimization of cost functions of the form 

V [Hy(t)llQ + II U (0IIr] (2.13) 

and 
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Q = diag[Q 0 ,QJ , Q > 0 


(2.14) 


provides a stable optimal closed-loop solution. 

3.0 SELECTION OF DYNAMIC SENSORS AND ACTUATORS 

3.1 Closed-Loop Input/Output Cost Analysis 

In order to write the expressions for the closed-loop input and output costs, it is first 
necessary to put the fully augmented system, under closed loop steady-state optimal 
state-estimate feedback control, in the following state space form: 

x(t) = Ax(t) + Dw(t) (3.1a) 

y(0 = Cx(t) (3.1b) 

v = £. V Q (t), V 0 (t) = y*(t)Qy(t) , (3.1c) 

where 

x T = [xW)> y T = Ly] T w T = [w T , v T ] (3. id) 

A BG DO CO Q ol \ W u] 

A= FM A+BG-FM ,D== 0 F ,C= 0G’^ = 0 1? ,W= jj* v (31e) 

G =-R~ l B T K, 0 = KA A t AT -KBR~ 1 B t K +C t QC (3.10 

F =[PM t +DU]V~\ 0 = [A-DUV~ X M]P +P[A-DUV- l M) r (3.1g) 

-pm t v~ x mp + dwd t - duv~ 1 u t d t 

For the system (3.1) the output costs V, y , defined by 
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V> = (l/2){£ a .0iy9y i )y i } 


(3.2a) 


are calculated as follows [2] 

V> = [CXC T QJ li 

where X is the steady state covariance satisfying 

0 = AX + XA t + DWD t 

and where the output costs satisfy the cost decomposition property 

iv>=v. 

i=l 


The input costs are defined by 


V i w = (l/2){E oo 0K (y /aw i )w i } 


and are found from [2] 

V ; w = [D T SDW] ii 

where S satisfies 

0 = A t S + SA + C t QC 

and where the input costs also satisfy the cost decomposition property 

n. 

IV, W = V. 


(3.2b) 


(3.2c) 


(3.2d) 


(3.3a) 


(3.3b) 


(3.3c) 


(3.3d) 


The input and output costs represent the in situ contributions that the noise inputs 
and the system outputs make in the cost function. We may also wish to know the amount 
by which the cost function will be reduced if a noise input is eliminated. This amount, 
AV*, is defined as 
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AV” = V - V Ri 


(3.4) 


where V Ri is the value of the cost function after the i^ noise input is eliminated, (but the 
controller is not redesigned) and AVj W is the cost reduction due to eliminating w ; . A 
positive value for AV i w indicates that elimination of the i^ input will reduce the cost, 
while negative AV; W indicates that a cost increase will follow noise elimination. It was 
shown in [1] that the AVj w may be positive or negative in the presence of noise 
correlation. Partitioning the matrices W and D facilitates direct solution for the cost 
reduction [2], yielding 

AVi W = 2V; W - d^Sd^ . (3.5) 


The closed-loop covariance X may be written 


X — - 


P+N N 
N N 


where P satisfies eqn (3.1g) and where N satisfies: 

0 = N(A +BG ) T + (A +BG )N + FVF T 


(3.6) 


(3.7) 


Also, S has the following form 


S = 


K+L -L 
-L L 


(3.8) 


where K satisfies eqn (3. If) and where L satisfies 

0 = L(A -FM ) 4 - (A -FM ) J L + G r RG (3.9) 


For notational convenience the steady state covariance X is partitioned as follows: 
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■ Xp x 12 x 13 

X=[P+N]= Xj 2 X. X23 (3.10) 

Xj3 x£ x s 

Using the notation of (3.10) and the special structure of the closed-loop system matrices 


in eqn (3.13) we write the following expressions for the output costs 

Vy^tCpXpC^Qplii i = 1, ny p (3.11a) 

V/ = [QXaCjQJji i = 1, • • • nu (3.1 lb) 

V i u = [GWG T R] ii i = 1, nu (3.11c) 

and for the input costs 

V i w = [D T (/i:+£-)^] ii i = 1, ■ • ■ nw (3.12a) 

V^ = [D T (K+L)DW] nw+inw+i i = 1, • • • nz (3.12b) 

V i Vo “ = (F T LFV] ii i=l, • nz (3.12c) 

and the input cost reductions 

AVj W = [D t (K+L)DW -D t LFU t ] ii i = 1, • • • nw (3.13a) 

AV i v " = [D J (K+L)DW -D t LFU t ] nw+inw+i i = 1, • • • nz (3.13b) 

&V?°“ = [F T LFV -F J LFV -F r LBU] d . i = 1, • • • nz (3.13c) 
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3.2 Dynamic Actuator Effectiveness Values 


Now that the closed-loop input and output costs have been determined for systems 
with dynamic sensors and actuators, it remains to use the CIOCA results to define 
expressions which reflect the effectiveness of each sensor and actuator in the cost 
function. This section defines the effectiveness values for dynamic actuators. The 
approach taken in [1] and [2] for non-dynamic actuators was to subtract the contribution 
the actuator’s noise in the cost function from the contribution of its control signal, and 
to label this difference the "effectiveness" of the i^, actuator, V* 01 . That is, 

V; act = V, u - AV* (3.14) 

This subtracts the "bad" from the "good" contributions of the actuator to measure its 
effectiveness. The results of applying (3.14) to sensor and actuator selection for a range 
of small and large scale examples in [2], [3], [4], [17] and [18] have demonstrated the 
utility of this approach. 

Extending the definition (3.14) for applicability to systems with dynamic actuators, 
we proceed as follows. In (3.1) there are two noise sources associated with each 
actuator: command noise, w u , which is filtered by the actuator dynamics; and output 
noise, w f , which is additive with the actuator output. Thus, the noise contribution 
associated with the i^ actuator is given by the sum of AVj W * and AVj Wf . 

The beneficial control cost for each actuator is not immediately evident. First, 
recall that it is the actuator output f(t), not its input u(t), which drives the system. Next, 
note that the contribution of the i^ actuator’s output in the cost function, V*, includes the 
effects of noise w ui . That is, even in the open loop (u = 0), V/ * 0 for [W u ] ;i > 0 with 
dynamics. Hence, to define the beneficial (control) portion of Vf it is necessary to 
subtract the portion of which is due to noise. This can not be accomplished exactly, 
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since the actuator command u(t) and the command noise w u (x) are correlated for t > x. 
An approximation is obtained, however, by solving for V* when u = 0 (that is, in the 
open loop). We define the contribution of w ui to V/ and the contribution of Uj to W( as 
follows, using the open loop covariance of the actuator states X a : 

[V i f ] w = [C a X a C a T QJ ii (3.15a) 

and 

[Vi f j u = V, f - [V ; f ] w = [C a (X.->yC jQJa (3-1 5b) 

where Xa solves 

0 = AaXa + XaAj + B a W u B a T . (3.1 5c) 

Finally, the input costs and the decomposition of the output cost Vf are combined in 
an effectiveness formula for dynamic actuators which is motivated by the results of [1] 
and [2]: 

V^ 1 = [V i f ] u - AV; Wf - AV*" . (3.16) 

Note that in the absence of command input noise, [V/] w and V*" are both zero. Also, in 
the absence of actuator dynamics, f ; (t) is equivalent to u ; (t). Thus the expression (3.16) 
reduces to the original effectiveness formula of [1] in the absence of actuator dynamics. 
Note also that (3.16) is applicable whether or not the actuator noise signals are correlated 
with other noise sources, and it is applicable to systems with actuator dynamics of 
arbitrary order. 


3.3 Dynamic Sensor Effectiveness Values 

Unlike the actuator noise, (which has a direct path to the output, independently of 
the conrollers influence) the noise associated with sensors reaches the system only 
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through the controller. Since the gains in the Kalman filter of the LQG controller 
represent an optimal trade-off of each sensor’s (beneficial) measurement information 
versus the (performance degrading) impact of its noise, then a AV ; V of large magnitude is 
indicative of a highly effective sensor. That is, the fact that a sensor’s noise is being 
allowed to heavily affect the cost means that its measurement information is even more 
critical to performance. For this reason, the following effectiveness formula for non- 
dynamic sensors, generalized to accommodate the possibility of noise correlation, was 
presented in [1]: 

V^ n =\AVf\. (3.17) 

For dynamic sensors there are two possible noise inputs associated with each sensor. 
As in the non-dynamic case, both noise inputs reach the system dynamics through the 
Kalman filter. Thus a straightforward extension of (3.17) to dynamic sensors is 

Vi sen = lAV^I + IAV i v “l . (3.18) 

Note that this formula is applicable in the presence of sensor dynamics of arbitrary order, 
and applies whether or not any of the noise sources are correlated with one another. 

This section concludes with the suggestion that (3.16) and (3.18) provide effective 
measures of the contribution of each actuator and sensor in a closed loop optimal LQG 
control (with sensor and actuator dynamics properly included). 


4.0 SELECTION OF DYNAMIC ACTUATORS FOR SCOLE 

In this section the actuator selection problem is solved for a model of NASA’s 
SCOLE (Spacecraft Control Laboratory Experiment) system. The SCOLE configuration 
consists of a flexible antenna suspended from the Space Shuttle cargo bay by a 130 ft. 
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flexible beam (see Figure 4.6). The effectiveness values for proof mass actuators 
(PMA’s) located along the beam are calculated and plotted versus position for both 
dynamic and non-dynamic actuators in order to evaluate the dynamic actuator selection 
method and to determine the effect of actuator dynamics on our results. 

4.1 PROBLEM DEFINITION 

A certain 2-dimensional SCOLE model includes four flexible modes and no rigid 
body modes [15-18]. Approximate open-loop mode shapes for the four flexible modes 
are presented in Figure 4.7, and the results of an open loop modal cost analysis are 
presented in Table 4.2. A detailed discussion of the model development is given in [16] 
and [18]. The two sensors retained in the model (using the CIOCA method of selection 
for non-dynamic sensors) are angular position and rate measurements located at the 
center of mass of the reflector [18]. Since there are no accelerometers presents, then the 
sensor and actuator noise is uncorrelated. Noise intensity data for the sensors is given in 
Table 4.3. 

The set of admissible actuators includes both a control moment gyro (CMG) located 
at the reflector center of mass and a set of PMAs distributed along the flexible beam. The 
actuator selection problem is to determine the optimal location for two PMA devices 
along the beam. To this end, the admissible set of PMAs was defined as 20 actuators 
spaced at distances of 6.25 feet apart on the 130-foot beam from a point 10.75 feet above 
the shuttle end of the beam to a point 129.5 feet from the shuttle. The PMA locations are 
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Table 4.2: 2-Dimensional SCOLH Elastic Modal Cost Analysis 


Mode# 

Frequency (Hz) 

Modal Cost 

Percent Total 

I 

.2&9E+00 

547E+01 

.466E+02 

2 

.164E+01 

626E+01 

533E+02 

3 

.497E+01 

.105E-01 

895E-01 

4 

124E+02 

. 1 K0E-04 

. 153E-03 

5 

.237E+02 

201E-06 

.171 E-05 

6 

.389E+02 

.692E-08 

.590E-07 

7 

580E+02 

470E-09 

.400E-08 

S 

81 0E *-02 

.SfffcE 10 

.43 i E 09 

Q 

. m^E *-03 

.74SE-1 1 

b'SE 10 

10 

. 1 39E+-03 

142E-11 

121E-K) 

11 

. 175E+03 

322E-12 

274E-11 

12 

215E-03 

H44E-P 

.723E 12 

13 

259E+03 

.243E-13 

.207E-12 


Table 4.3: Noise Specifications for SCOLE Actuators and Sensors 

Actuators 

Dynamic Noise Noi^e 

T>pe Range !nten,:ly Type I..:crts:i 


PM a 10 lb 0001 <Ib) 2 Accelerometers v # ■= .0025 (deg/sec*) 2 

CMG 10 3 ft-lb 10,000. i. ft— lb) 2 Angular Pos t:on I O' * i deg) 2 

004 Meg/sec) 2 
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thus selected by evaluating the relative effectiveness of each of the 20 PMA locations. 


4.2 RESULTS FOR NON-DYNAMIC ACTUATORS 

The PMA selection problem for non-dynamic actuators was solved first, for later 
comparison with the dynamic actuator selection results. In all cases (dynamic and non- 
dynamic) the actuator effectiveness values are calculated following controller design 
which achieves a specified output variance and minimizes the amount by which the 
actuators exceed their specified variances. This type of controller is designed by an 
iterative selection of the control and output weights using the Output Variance 
Assignment (OVA) algorithm (DeLorenzo and Skelton, [3]). The variance specification 
for each actuator was equal to 10 times the intensity of its noise (see Table 4.3). 

The actuator effectiveness values based upon standard Closed Loop Input/Output 
Cost Analysis (CIOCA) [2] for non-dynamic actuators are presented in Figure 4.8. The 
figure portrays PMA effectiveness results for four different controllers, each achieving a 
different steady-state line-of-sight (LOS) error variance. The results provide a vivid 
illustration of how the controller objectives can profoundly influence the actuator 
selection results. For lower gain controllers (lower LOS error) the theory determines that 
the upper tip is the most desirable PMA location. However, as the gain increases 
(controller designed for smaller LOS error) the center of the beam becomes the optimal 
location. 

The results of Figure 4.8 are readily explained via modal analysis. The mode shape 
figures for the four flexible modes retained in the 2-dimensional SCOLE model were 
presented in Figure 4.7. Recall that mode #1, which accounts for 46.6 percent of the 
open loop modal cost, has a maximum amplitude at the reflector-end tip of the beam (i.e., 
at 130 ft.). Mode number #2, which accounts for 53.3 percent of the open loop modal 
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cost, has a maximum amplitude near the 90 ft. point. And mode #3, which accounts for 
only approximately 0. 1 percent of the open loop modal cost, has a peak amplitude near 
the center of the beam. 

Next note from Figure 4.8 that as the gain is increased in order to achieve a smaller 
steady state LOS error variance, the most effective location for PMAs shifts from the tip 
of the beam to the midpoint. This corresponds to a shift from the peak of mode #1 to the 
peak of mode #3. The shift occurs even though with higher gain the noise in the PMAs 
near the beam midpoint becomes the most detrimental to performance (Figure 4.9). In 
fact, Figure 4.9 indicates the reason for the shift in optimal PMA location: with higher 
gain the third mode becomes the least damped by the control of the CMG, and becomes 
therefore a significant mode to be controlled by the PMAs. Figure 4.10 and Table 4.4 
indicate the motion of the closed-loop eigenvalues from their open loop locations under 
varying levels of gain (output performance). 

Since the control cost of each PMA (V“ = E^qu 2 ) is equal to its effectiveness value 
V® ct minus the cost contribution of its noise, V*', then it is clear from Figures 4.8 and 4.9 
that the PMAs are being used primarily to control mode #3 (i.e., near the middle of the 
beam). However mode #3 is the most lightly damped mode in the closed loop. This is 
true in spite of the fact that in all cases the input variances of both the CMG and PMAs, 
when normalized by their variance specifications, are of like order of magnitude (see 
Figure 4.1 1). 

The results demonstrate the interesting result that while the PMAs are being used at 
a level similar to the CMGs (in relation to their specified variance levels), they 
nontheless make a small contribution to the closed-loop modal damping. This claim is 
verified by deleting all PMAs from the system and again using OVA to achieve a 
specified LOS error of 0.1 (arc sec) 2 , and comparing the resulting closed-loop modal and 
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Tabic 4 4: Modal Characteristics of l A+BGj as a Function of Output Performance 
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Loop 

.1015 

(deg) 2 
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(deg) 2 
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(deg) 2 
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1.8 

3.14 

3.14 

3 14 

3.14 

3.135 
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.005 

,028 

,027 

.0255 

.025 

.01 

Tj (sec) 

mi 
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12 5 

12.74 

31.9 
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10.3 
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160. 
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.015 

009 

023 

o^fr/s) 

31 
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28.8 

28.8 
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.005 

.0217 
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.009 
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3.16 
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performance data with that obtained from a full set of PMAs and an output variance of 

0.1 (arc sec) 2 (see Table 4.4). 

4.3 RESULTS FOR DYNAMIC ACTUATORS 

In this section we add actuator dynamics to the SCOLE model and then re-solve the 
actuator selection problem solved above. The actuator dynamics are given in NASA’s 
original SCOLE document [15] to be first-order with a time constant of 0.1 seconds. 
That is, for each actuator (both PMA and CMG) the response of the actuator f,(t) to its 
input signal u,(t) is governed by 

f^sVu^s) = [l/(.ls+l)] (4.7a) 

or 

= — 1 Of i + 10u ; . (4.7b) 

There are several possibilities for the characteristics of the white noise associated 
with the actuators; white noise may be an input to (and thus be filtered by) the actuator, 
or it may be additive with the actuator output (thus unfiltered), or both. In this example 
four different actuator noise models are considered. Recalling that the non-dynamic 
actuators had additive white noise with intensity W, the following noise cases were 
studied for dynamic actuators: 

1. white actuator input noise of intensity W u = W; 

2. white actuator output noise of intensity W f = W; 

3. both input and output noise, each white and of intensity W; 
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4. both input and output noise, each white and of intensity W/2. 

The sensors are assumed non-dynamic (without phase lag). 

First we examine the effect of actuator dynamics on the maximal theoretically 
achievable accuracy. From [3], the lower bound y * on the steady-state variance of the i^ 
output is given by 

y* = [CPC T ] H i = 1, .... ny . (4.8) 

The values of the lower bound on the LOS error for the fourth-order 2-D SCOLE model 
under study were calculated for the four different actuator noise cases listed above, as 
well as for the non-dynamic actuator model examined earlier. The results are shown 
below. 


Table 4.5: Maximal Accuracy for Different Actuator Noise Cases 


Noise Case 

No Dynamics 

W U =W 

WpW 

Wf=W u =W 

Wf=W u =W/2 

Max. Acc. 

.086921 

.0691 

.086921 

.10072 

.07926 

(arc sec) 2 







From Table 4.5 it is clear that the addition of actuator dynamics along with 
retention of the white noise input to the system states (actuator output noise only, Wj=W) 
does not change the theoretical maximal accuracy; that is, y* is equal for the non- 
dynamic and the W f = W case. Also from the table, filtering of the actuator noise by 
passing it through finite actuator dynamics clearly improves the maximal accuracy. 
Finally it is noted that for case (3), W f = W u = W, the minimal LOS error is greater than 
that obtained by all but one of the controllers in the non-dynamic case. Thus for 
purposes of comparison only cases (1), (2) and (4) are studied in further detail. 
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For each of the three actuator noise cases a controller was designed (using OVA) 
which assigned the steady state LOS error variance to 0.1 (arc sec) and minimized the 
sum of the normalized actuator variances among those actuators whose variances exceed 
their specifications (normalized variances greater than unity). For each final controller, 
the dynamic actuator effectiveness values for the PMAs are plotted in Figure 4.12 versus 
the actuators’ position along the 130 ft. flexible beam. 

For each of the noise cases the most effective actuator location is toward the 
reflector-end of the beam, with the highest effectiveness values corresponding to 
actuators located at the beam tip. Recalling the mode shape figures for the open loop 
flexible modes, the results in Figure 4.12 indicate that the PMAs are used by the optimal 
controller primarily for control of mode #1, which accounted for 46.6 percent of the open 
loop modal cost. It is interesting to compare Figure 4.12 with the plot of effectiveness 
values for non-dynamic actuators (Figure 4.8); note that the most effective non-dynamic 
actuators for the controller which achieved LOS error = 0.1 (arc sec) 2 were lo:ated near 
the center of the beam (70 ft from the shuttle). Hence, the optimal beam locations for 
PMAs in controllers which are achieving the same output performance are affected by 
the actuator dynamics. 


CONCLUSIONS 

The Closed-Loop Input/Output Cost Analysis (CIOCA) method of sensor and 
actuator selection (SAS) has been extended for application to systems with dynamic 
sensors and actuators — that is, systems in which the response of the sensors and 
actuators to their inputs is not instantaneous but governed by deterministic dynamics. 
The extended SAS method is applicable to systems in which the deterministic sensor and 
actuator dynamics are of arbitrary order. Application to simple numerical examples 
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demonstrates the utility of the SAS method. The examples also demonstrated that even 
uniform sensor dynamics can affect the optimal selection of sensors. Application of the 
actuator selection method in detail to NASA’s SCOLE space structure demonstrated that 
even uniform actuator dynamics can affect the optimal selection of actuators. 
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Proof of Theorem 1 


Let a = (nx p +nx a ). The composite system (2.7e) has a controllability matrix 

W c 6 R a *( a ’ nu ) of the following form 

0 B 0 C a B a (ApB 0 C a B a +B 0 C a A a B a) 

B a A a B a A a B a 

• • • (A p a - 2 B 0 C a B a +A p a - 3 B 0 C a A a B a + • • • +B 0 C a A«- 2 B a ) 

(A1) 

Now noting that C a A a B a = Mj is the ith Markov parameter for the system of actuator 
dynamics (2.7c,d), W cl may be rewritten 

W cl = tO B 0 M o (ApB^+B^j) (A 2 B 0 M 0 +A p B 0 M 1+ B 0 M 2 ) • • • 

• • • (A p a - 2 B o M o +A p a - 3 B 0 M 1+ • • • +B 0 M a _ 2 )] (A.2) 

The columns of W c span the controllable subspace of the composite system. Linear 
independence of all the rows in W c implies full controllability of the composite system. 
However, controllability of the original s\ stem states, x p , requires only that the columns 
of W c l span the state space for x 0 . This in turn will hold if and only if the matrix W cl has 
rank nx p . 

The proof of the Theorem begins with the proof that (A. 3) implies (A. 4): 

(det M k * 0, Mj = 0, i = 0, 1, .... k-1 } (A.3) 

rank[W cl ] = n x ^ (or range space of W c[ has dimension n x ) (A. 4) 

Note that the last block of W el has the form 

w cl(0) = (A?- 2 B 0 M 0 +A“- 3 B 0 M 1 + • • • + B 0 M o _ 2 ) (A. 5) 

Now let k ^ a-2 be the index of the first nonzero Markov parameter, M k . (In this case 


W c = 


w, 

w 


cl 


c2 
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the first k+1 blocks of W cl are zero.) Next, use is made of two results from linear 
algebra ( "R [K]" denotes "range space of K"), 

{detK*0} => [R[JK)=R[J]} (A.6a) 

R[J+K]czR[J] + R[K] (A. 7b) 

(where "c" means "is contained in") to demonstrate the following results which hold 
when M k is nonsingular 

R [B 0 ] = R [B 0 M k ] (A. 8) 

R[B o A p B o ]=R[B 0 ]+R[A p B o ] 

~ R [B 0 M k ] + R [A p B 0 M k ] 

= R [B 0 M k ] + R [A p B 0 M k + B 0 M k+1 -B 0 M k+1 1 
<= R [ 0 oM k ] + R [A p B 0 M k + B 0 M k+1 ] + R [B 0 M k+1 ] 

= R [B 0 M k ] + R [A p B 0 M k + B 0 M k+1 ] 

R [B 0 A p B 0 ] c R [B 0 M k A p B 0 M k + B 0 M k+1 ] (A.9) 

Eqns (A. 8) and (A.9) lead by induction to the main result 

R [B 0 ApB 0 ••• A “ -k_2 B 0 ] c R LB 0 M k A p B 0 M k + B 0 M k+1 ••• 

A p a - k - 2 B 0 M k + • • • B 0 M a _ 2 ] 

that is, 

ft |B 0 A p B 0 ■ ■ ■ A“- k ' ! B„] c ft [W c , 1 . (A. 10) 

254 


I 



Condition (A. 3) leads to (A. 10). Thus, given (A.3) together with (A p ,B 0 ) 
controllable, the columns of W cl are guaranteed to span the nxp-dimensional state space 
for Xp as long as 

a-k-2 > nx Q -l . 

that is, as long as 

k < nx a - 1 . (A.l 1) 

In fact, the index k of the first nonzero Markov parameter for the system (2.7) will 
always satisfy (A. 11). To show this, simply note that by observability of (A a ,B a ), the 
observability matrix W oa for (2.7) has full column rank: 

rank(W oa ) = nx a (A. 1 2) 

From (A. 12), 

{W oa B a = 0} => (B a = 0) => {Contradiction of (A a ,B a ) controllable} (A. 13) 

Thus, 

W oa B a = [M g t , Mj , .... M 7 ^] 7 * 0 (A.14) 

and so the validity of (A.l 1) is guaranteed for (2.7) completing the proof that 

{ [(A p ,B 0 ) controllable] & [M k * 0, M* = 0, i = 0, 1 , ..., k-1 ] } (A. 1 5) 

=* {x p controllable u} 

The usefulness of (A. 15) stems from the fact that by meeting certain mildly 
restrictive conditions the actuator dynamics can be guaranteed not to destroy 
controllability of the original system states x p , regardless of the pole/zero location for the 
plant. 
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It remains to prove the equivalence of the condition (A.3) and the requirements on 
the individual actuators’ numbers of poles and zeros. First, note that since each actuator 
is a single input, single output (SISO) system, then the Markov parameters for the 
lumped actuator dynamics (2.7) are diagonal matrices of the following form: 

M; = diag [m H , m 2i , m 3i , m^J (A. 16) 

where m^ is the (scalar) i^ Markov parameter for the actuator. Thus the condition 
(A.3) is met if and only if the index i of the first nonzero markov parameter is equal 
among all the actuators. 

The input/output transfer function for any n^, order SISO system has the form: 

T(s) = (c n _ l s n - 1 +c n _ 2 s n " 2 + • • • + c 0 )/(s n +d tl _ 1 s n-1 + • ■ • +d 0 ) (A. 17) 

The scalar Markov parameters m { for the SISO system with transfer function (A. 20) may 
be shown to be given by: 


m l = c n-2-. d n-l n o 


n 2 — c n-3 — ^n-2 no^-l 0 ! 

n n = c o~di -d2 n l ~ -dn-iV-2 


(A. 18) 


From (A. 18), n t is the first nonzero Markov parameter for a system when the 
number of zeros in its transfer function is 

z = n - i - 1 (A. 19) 

Letting nj and Zj equal the number of poles and zeros for the actuator, respectively, 
(A. 19) yields the conclusion that 
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{det M k ^ 0, Mj, = 0, i = 0, 1, k-1 ) ** {(nj-Zj) = (rij— Zj) V i,j,e (l,2,...,nu)}(A.20) 


Thus it is concluded that 

{(n j -z j ) = (n r z i )V i,j,e(l,2,...,nu)} => {x p is controllable u}. 


(A. 21) 
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